from math import *
x=[24.533,16.033,35.95,7.85,0.583,16.867,6.95,0.333,18.25,10.583,14.8,7.333]
y=[0.34,0.264,0.443,0.19,0.126,0.271,0.182,0.123,0.283,0.215,0.252,0.185]
n=len(x)
vmx=sum(x)
vcx=sum([w*w for w in x])
vmy=sum(y)
vcy=sum([w*w for w in y])
vcx=vcx-vmx**2/n
vcy=vcy-vmy**2/n
xy=sum([w*z for w,z in zip(x,y)])-vmx*vmy/n
b=xy/vcx
a=(vmy-vmx*b)/n
d=b*xy/vcy
print("régression lineaire y = a + b * x".center(40))
print("{:s}{:>10.3g}{:s}{:>10.3g}{:s}{:>10.3g}{:s}".format("a : ",
    a,"\nb : ",b,"\nd : ",d,"\n\n      x       y-calc   y-donné"))
for i in range(0,len(x)):
    
    print("{:>10.3f}{:>10.3f}{:>10.3f}".format(x[i],a+b*x[i],y[i]))
#print()

